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Abstract 

At the example of two coupled waveguides we construct a periodic second order dif- 
ferential operator acting in a Euclidean domain and having spectral gaps whose edges arc 
attained strictly inside the Brillouin zone. The waveguides are modeled by the Laplacian in 
two infinite strips of different width that have a common interior boundary. On this com- 
mon boundary we impose the Neumann boundary condition but cut out a periodic system 
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of small holes, while on the remaining exterior boundary we impose the Dirichlet boundary 

condition. It is shown that, by varying the widths of the strips and the distance between 

the holes, one can control the location of the extrema of the band functions as well as the 

number of the open gaps. We calculate the leading terms in the asymptotics for the gap 

lengths and the location of the extrema. 
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Introduction 



It is well known that a large class of periodic operators have a band spectrum, i.e. the spectrum 
is the union of finite segments. One of the important questions in the study of such operators 
concerns the existence of spectral gaps, i.e. open intervals lying outside the spectrum but whose 
edges are in the spectrum. The study of this property is motivated by various applications, in 
particular, by the problems arising in the theory of photonic crystals fjQ. 

The one-dimensional Schrodinger operator with a periodic potential (Hill operator) is a clas- 
sical example of a periodic operator. The only possible gap edges for these operators are the 

CN ' values of the energy for which the eigenvalue equation has a non-trivial periodic or anti-periodic 

solution. While this reflects some specific features of the operator, the same was implicitly 
assumed to be true for more general operators; in those cases the periodic and anti-periodic 

S^ | eigenvalues are replaced by the values of the band functions at some special values of the quasi- 

momentum (see below). A discussion of analytic aspects of this property appeared just quite 
recently in the works [21[3], which served the motivation for the present paper. The work [3] 
contains, in particular, several examples of quite sophisticated multi-dimensional combinatorial 
operators whose gap edges are attained at intermediate values of the quasimomentum, and the 
work [2] is devoted to a more detailed analysis of Z-periodic operators. 

It should be emphasized that similar effects, like attaining extremal values at intermediate 
values of quasimomenta and the presence of multiple extrema, are known in the physics literature 
for a long time as split band edge, and these effects play an essential role, in particular, in slowing 
waves in photonic crystals, see e.g. [IHS]. At the physics level of rigor it was noted that the split 
band edges can appear when coupling several Z-periodic structures with different characteristics, 
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cf. [6]. In the present work we give this assertion a rigorous sense at the example of two non- 
symmetric waveguides. We consider the systems consisting of two two-dimensional strips of 
different width coupled by a periodic system of windows. The Dirichlet boundary conditions are 
imposed at the exterior boundary, and the Neumann boundary conditions at the interior one. 
It is shown that, if some inequalities between the period and the waveguide widths hold, under 
the assumption that the windows are small, the associated operator (Laplacian) has gaps whose 
edges are attained at intermediate values of quasimomentum. We note that we use explicitly 
the Z-periodicity (i.e. periodicity in one direction and compactness in the others), and we do 
not expect the results to be valid for systems with multi-dimensional periodicity (see e.g. the 
review [7] for possible ways of gap opening in multidimensional systems). 

The questions of gap opening in coupled and perturbed waveguides (in situations which 
are different from the one under consideration) were already studied in numerous works. In 
particular, the existence of gaps in the spectrum of a periodically curved waveguide was shown 
in [8], and the work [9] deals with the study of the asymptotics of the bands and the gaps 
of curved waveguides with respect to the cross-section size. Several works studied by various 
methods the waveguides represented as a chain of domains coupled by windows or links |10H13j ; 
in that case the number of gap grows to infinity. The works [14H16J are the closest to us 
from the point of view of methods, they show the gap opening in cylindrical waveguides with 
periodic localized perturbation at the boundary. Furthermore, there were some works in the 
physics literature |17H19j . dealing with the study of gaps in waveguides coupled or perturbed by 
periodically varying boundary conditions. The main contribution of the present work consists 
in a rigorous description of the split band edge effect in coupled waveguides, which, to our 
knowledge, was not found earlier. 

The paper is organized as follows. In section [1] the rigorous problem setting is given and the 
most important result, Theorem ll.il is formulated, whose proof is given in the two subsequent 
sections. In section [2] we are concerned with an upper bound of the shift of band functions of 
uncoupled and coupled waveguides. In section[3]the crossings of the unperturbed band functions 
are analyzed in greater details, and the asymptotics of the perturbed band functions near these 
crossings is given. 

1 Problem setting 

For the sake of convenience let us denote <i+ := n, d~ = d > 0, and, furthermore, let h and e be 
parameters satisfying ^ e < h. In R 2 we consider infinite straight strips 

11+ :=Rx (0,d+), n_ :=Rx (-d_,0), II := R x (-d_, d+), 

and introduce the operator H e acting in 1/2(11) as the Laplacian with the Dirichlet boundary 
conditions at the lines X2 = ±d± and the Neumann boundary conditions at the line X2 = 
except at the intervals w^ := (2nh — e, 2nh + e) x {0}, n£2. It is clear that for e = we just 
have the equality H° = H + © H— , where H± are the Laplacians in the waveguides TT-t- with the 
Dirichlet boundary conditions at the lines %2 = ±d± and the Neumann boundary conditions at 
the line X2 = 0. 

Therefore, the operator H° is the Hamiltonian of the system consisting of two non-interacting 
waveguides 11+ and n_, while H £ describes the same waveguides coupled by the 2/i-periodic 
system of windows u^ of length 2e, see Figure [TJ 

The present works analyses the situation of non-symmetric waveguides: d ^ it. Without loss 
of generality we assume 

d<TT. (1.1) 
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Figure 1: The waveguides Il + =Rx(0, n) and II_ =Rx (— d, 0) coupled by a periodic system of 
windows. The parts of the boundary on which the Dirichlet and Neumann boundary conditions 
are imposed are marked by the symbols D and N respectively. The elementary cell H h is 
highlighted. We denote Il£ := U h n II±. 



The spectra of the operators H± are well studied. There holds 
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hence, under the assumption (jl.lj) . the spectrum of H° is the half-axis [1/4, +oo). Moreover, 
the spectra of H± and H° are absolutely continuous. 

Let us denote by II the elementary cell given by II := II fl {x = {x\, X2) '■ \x±\ < h}. The 
periodicity with respect to x\ and the application of the Bloch-Floquet theory allows one to 
decompose H £ into a direct integral: 



; / H £ (k)dk, 

' (— 7T,7r] 



where H e (k) is the operator acting in L2(Ii h ) by 
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with the boundary conditions 
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where 



r D := {x = (xi, d+) : \x\\ < h\ U {x = (x\, — gL) : |»i| < /i}, 
T N ' e := {x = (xi, 0) : £ < \x\\ < h}. 
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The variable k in the preceding formulas is usually referred to as the quasimomentum or the 
Bloch parameter. 

The spectrum of each of the operators H e (k) is discrete. Let us denote by Ef(k), I G 
IN, its eigenvalues taken with their multiplicities and enumerated in the non-decreasing order, 
Ef(k) ^ Ef +l {k). The functions k \— > Ef(k), called band functions (or dispersion laws, or 
dispersion relations), are 27r-periodic, continuous and even (in virtue of the real-valuedness of 
the Laplacian), and the spectrum of H e is nothing but the union of the ranges of these functions, 

a(H*)=\jBf, Bf = Ef{(-^]). 

The segment Bf =: [a7,aii~] is called Ith. band of the operator H e . If for some I G IN one has 
af < ctf+i' then the interval (a^~, a^ +l ) is called a gap of the operator H e . 

The central question of our work is the question on the existence of gaps in the spectrum 
of H e for small e, as well as the question on the values of the quasimomentum k at which the 
extremal values a ; of the dispersion laws are attained. In particular, the attaining the extremal 
values outside the particular values k = and k = ±7r will be put into evidence. 

To formulate the main results we need some additional notation. Let us denote Ilij_ := 
n± n {x : \x\\ < h}. In the spaces ^(II^.) consider the operators H±(k) given by the differ- 
ential expression (|1.3p with the boundary conditions (|1.4p . where in the Neumann condition 
(|1.4bp we formally assume T N ' := {x = (xi,0) : \xi\ < h}. Analogously to (|1.2|) . one has the 
decomposition 

H±~ H±{k)dk. 

The operators H±(k) have discrete spectra, and their eigenvalues and eigenfunctions are 
calculated explicitly. The eigenvalues are 

, „. /A; + 27rm\2 ( it \ 2 ( 1\ 2 , -, , N 

E±, p {k) = { 2h ) +(^j (p+-) , meZ, P GZ + ^1NU{0}, (1.5) 

for k ^ {0,7r} the eigenvalue E^ (k) is simple with the eigenfunction 

while for k G {0, ir} the respective eigenvalue is double. The spectrum of H°(k) is exactly the 
union, taking into account the multiplicity, of the eigenvalues E+ [k) and i?~ p (k) for all m G Z, 

Let us note that the multiplicity grows infinitely at high energies: each point of the spectrum 

E of the operator H° belongs the the ranges of k(E) function E^ , and k(E) > +oo. Our 

main result is contained in the following theorem and describes the behavior of the spectrum 

of H e in the neighborhood of the points E with k(E) ^ 4. Note that the inequality k(E) ^ 4 

9 
implies automatically E < — . 

Theorem 1.1. Assume that the inequality II 1.1]) holds and that there exists a point ko G (0,vr) 
and m,n G Z for which one has 

E := E+ fi (k ) = E mfi (k ) < j, (1.7a) 
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Figure 2: A schematic view of two possible configurations of band functions for H°. The solid 
lines correspond to the functions E+ , the dashed ones correspond to the functions E+ p . (a) 
The intersections of the band functions satisfying the assumption of Theorem 11.11 are marked 
by fat points, (b) There are no points satisfying the assumptions of Theorem ll.il 
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_J£ (W . _££(„) < o. (i.7b) 

Then, for sufficiently small e, there is a gap (a/(e),a r (e)) in the spectrum of the operator H £ 
whose endpoints have the asymptotics 

aj{£) = Eo ~ih[kF\ +0{lln£r2) > i = / ' r ' (L8) 

where 



n -.= 2VC(i -P 2 ) - /3(C - i) - C - i, 

r r := -2 VC(1 - P 2 ) -/?(C-1)-C-1, (i.g) 

_ 7r _ 7r(m + n) + k 

<i' 7r(n — m) 

The band functions Ef(k) and E £ (k) describing the points of the spectrum to the left, respectively 
to the right, from the gap (a;(e), a r (e)) ; attain their extremal values ai(e) = maxm.^ Ef(k) and 
a r (e) = min^^ E^(k) at the points ±ki(e) and ±k r (e), respectively: 



Ef(±k(e)) = ai(e), E e r { ± k r (e)) = a r (e) 



One has the asymptotics 



k l/r (e) = k + ^ + O(\lne\- 3 / 2 ), 
' me 

ai vr(n - m) V 1 - /3 2 2vr(n - m) ' (1-10) 



/9h C (l-C)/i 

o"r = — ; r \ / ^ ^7 + 



ir(n — m) y 1 — /3 2 27r(n — m) 

Let us remark that one can use the methods of the works |2U[ Ch. 9] , J2ip23| and proceed with 
a deeper analysis in order to obtain further terms of the asymptotics (jl.lOp as some powers of 
1/ In e. Moreover, one can show some analyticity properties with respect to 1/ In e, see [2U|J22(I23| . 

Note that the theorem stated does not guarantee that the inequalities (|1.7p hold for at least 

one combination (m, n, ko). Moreover, one can easily choose the parameters in such a way 

that this condition will not be satisfied at any point; for example, this can be obtained by 

taking h sufficiently small (see Fig. [2] for illustration of possible configurations). On the other 

hand, to satisfy the conditions (jl.7p in at least one point it is sufficient to have the inequality 

9 
Eq (0) < E_ 10 (0) < -, see Fig. [2jb). Rewriting this inequality with the explicit expressions 

(11.51) we arrive at 



Corollary 1.2. Let the conditions 

h > —= and — , = < a < ir 

hold. Then there exists ko E (0, it) for which |j.7p holds with n = — 1 and m = 0, and the operator 
H £ , for e sufficiently small, has a gap lying in an 0(ln~ e) -neighborhood of the respective 
point Eq, and the adjacent band functions attain the respective extremal values in an C(ln~ e)- 
neighborhood of k^ . 
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Figure 3: Gap opening in the situation of Theorem 11.1 



Figure [3] illustrates the behavior of the band functions near the points at which the assump- 
tions of Theorem 11.11 Note that even though the technique of proof of Theorem 11.11 can be 
transferred without significant additional efforts to the case k G {0, tt} and allows one to clarify 
the behavior of the band functions at such crossings (domain f2 on Fig. [3]) by some additional 
computations, these parts are overlapped by the other parts of the graphs when projecting to 
the axis E, hence these crossings do not participate in the gap opening. 

We remark that by a suitable choice of the parameters d and h one can satisfy the assumptions 
of Theorem II. II at any prescribed value of quasimomentum Uq G (0, tt). In a similar way one can 
obtain any prescribed number of crossings in the interval (1/4, 9/4), and each of them will satisfy 
the conditions (|1.7p . Each of these crossings will generate its own gap by Theorem 11.11 Hence 
by fitting the parameters d and h one can open any prescribed number of gaps on the interval 
(1/4,9/4), and each of them is covered by Theorem 11.11 The respective values fco associated 
with different gaps will, in general, also be different. 



2 Estimating the shift of band functions 

One of the main ingredients of the proof of Theorem 11.11 comes from a certain rather general 
estimate of the difference \Ef(k) — Ef(k)\. This is discussed in the present section. The estimate 
needed is contained in the following Lemma. 

Lemma 2.1. For any I G M there exists C\ > such that, for e > sufficiently small, and all 
k G (—71", n] one has the inequality 



0^Ef(k)-E?(k)^ 



Q 
In el 



(2.1) 



For the proof of the required estimates, like how it was done in similar situations in [lOl lllj. 
we use the max-min principle. Let us consider the quadratic forms a® and a £ k associated with 
the operators H°(k) and H £ (k) respectively. The standard constructions show that the form a® 
is given by the expression 



a° k (u,u) 
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X « 7T7- J^l,^, 

ox\ 2h; 
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du(xi,x 2 ) 



dx 2 



dxi dx2 



on the domain D° consisting of the functions u G W 2 l (n+) © W^iJ^l.) satisfying the boundary 
conditions 

u\ rh = 0, V,\ Xl =-.h = U\xi=h- ( 2 - 2 ) 



D 



In turn, the form at is the restriction of a^ to the set D 6 consisting of the functions u G D° 
satisfying, in addition to (J2.2J) . the matching condition on the window 

^ + Ug = -""Ug, (2.3) 

where u stands for the restriction of u onto 11^ . 

According to the max-min principle |24} Th. XIII. 2], for the eigenvalues Ef(k) one has the 
equality 

Ef(k) = max min P$^ , I G IN, e > 0, (2.4) 

L e s ; _ lue ^nL|M|2 2(nh) 

where we denote by S m the set of all closed subspaces of L 2 (H h ) of codimension m, (i.e. dim L 1 - = 
m for all L G S m ). As D £ C D , the left-hand part of the inequality (|2.ip follows directly from 

dZ3D- 

To prove the right-hand part of (|2.1[) we use the approach similar to the one used in [10|]11|. 
Introduce first a family of cut-off functions. Let ip : [0, +oo) — > R be an infinitely smooth 
function with i/j(s) = 1 for s ^ ^ and ip(s) = for s ^ 1. For e > let us introduce the 
functions $ £ : R 2 — > R by the equalities 



$ £ (x) = v( 



In \x\ 



lne 



x + 0, * e (0) = 0. 



Note that for any function v G D° the product u = <& £ v for e sufficiently small is still in D° and, 
moreover, vanishes near the window ojq. Therefore one has the equality ()2.3p and the inclusion 
u G D e . Moreover, denoting B\ = max sS R |t//(s)|, for e sufficiently small we have the estimates 



V$p \dx 



u h 



2vr 
In el 



v^ 



V 



, flnr 

lne 



dr < 



2TiBiy/e 
I In el 



(2.5) 
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(2.6) 



Now let I G IN. Let us pick the eigenfunctions vh, j = I,..., I, of the operator H°(k) 
associated respectively with the eigenvalues E^{k), i.e. H a (k)vh = E^{k)vh, in such a way that 
the normalization conditions 

I k k\ I ' 3 = ^' 

\ u ji u slL 2 m h ) = \ „ . , 

[0, Jy^S. 

hold. Note that the restrictions of the eigenfunctions u^ onto the rectangles 11^. are either zero 



functions or given by the explicit expressions (|1.6j) . It follows that one can choose a constant 
i?2 > such that 

sup \Uj(x)\ + \VUj(x)\ <B 2 , j = l,...,l, fe€(-7T,7r], 
which can be used to prove that there exists a constant S3 > for which 



($ £ 4,$ £ ^) L 2 (nh ) " («J,«J>£3(nfc) 



<B 3 8 



(2.7) 



for all j, s = 1, . . . ,1 k G (— 7r,7r]. In particular, it follows from this last estimate that the I 
functions $ e iC , j = 1, . . . , I are linearly independent for small e > 0. Moreover, as noted above, 
each of these functions belongs to D £ . 



Let us consider now an arbitrary subspace L E Si-±. From the dimension considerations, 
the intersection D £ n L should contain at least one linear combination v of the form 

I 
« = ^6 i * e «}, b = (b J )e€ l \{0}, 

and, due to (|2.7p . for any vector b there holds 

with some constant B4 > which is independent of e and k. By elementary computations we 
obtain 
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it follows, by using (|2.5p and (|2.6p that for all j, s = 1, . . . , I and k G (— 7r, 7r] there holds 
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with some positive constants B$,Bq,B7 independent of k and e. 

By noting that a^{u k ,u k ) = E^{k)(u k ,u k ) £,am'») an d combining all the above estimates we 
arrive at 



a° fc (^K5>^)|6/ + 
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where the constant B$ > does not depend on 6, e and fc. Now, due to Ef(k) 
maxjgji ^ Ej(k), we have 
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with some constant C\ independent of b,e,k. As the subspace L G S/„i was arbitrary, the 
right-hand side of ()2.ip follows from the max-min principle (|2.4p . 



3 Analysis near the crossings of band functions 

In this section we proceed with a more detailed analysis of the eigenvalues Ef(k) near k = ko 
and conclude the proof of Theorem 11.11 

It follows directly from the definition of the functions E m that if the condition (|1.7|) holds 
for some m, n, ko, then it is also valid for —m, —n, —ko, and the value Eq remains unchanged. 
Furthermore, there are no other points at which (|1.7ap holds with the same value Eq. Therefore, 
it follows from Lemma 12.11 that for any A > and sufficiently small e there exists a constant 
to = to(A) > for which the inequality 

dist( a(H J k)),E ) > n r for \k - k \ ^ -p- and \k + U > 7^-7 (3.1) 

|lne| |l n£ | U ne l 

holds. Therefore, for those k which are sufficiently far away from ±fco, the spectrum of H £ {k) is 
separated from the point Eq by a distance of at least 0(\ lne| -1 ). For the proof of theorem it is 
now necessary to investigate the behavior of the spectrum of the operator H e (k) as \k — ko\ ^ 
to\ lne| _1 and \k + ko\ ^ to\ hie| _1 . Such a study will be done in a neighborhood of the point 
ko only. For the other case, the results can be easily obtained using the obvious fact that the 
replacement k h-> — k in the operator H £ (k) does not change the eigenvalues, and the associated 
eigenfunctions are complex conjugate to the original ones. 
We will consider the values k having the form 

t 
k £ = k + - — , 
me 

where t is a real parameter running through a segment [— to, to]- By Lemma 12.11 the operator 
H e (k £ ) has exactly two eigenvalues converging to Eq for e — > +0. Denote them by E e (t), 
j = 1, 2. We will assume the ordering E £ (t) ^ E £ (t). 

To construct the asymptotics of E £ (t) for e — > we will use the methods of matched 
asymptotic expansions [25]. This will be done in two steps: first we construct the formal 
asymptotic expansions, and then they will be rigorously justified. 

3.1 Constructing formal asymptotics 

We will construct the asymptotics of the above eigenvalues in the following form: 

E^(t) = E + ^f + ..., (3.2) 



where /i"'(t) are some constants and should be determined. 

wi: 
of H £ {k £ ) associated to the eigenvalues E £ will be sought in the form 



(i) 
Outside a small neighborhood of the window Wq the asymptotics of the eigenfunctions ips 
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with 



UiUfttiix), * 2 >0, ^(-)=e 1 ^^-cos y , 

\aPfi{x), x 2 <0, ^- (x) = cos !g 5 
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where a± are constants which cannot vanish simultaneously. Everywhere below we call the 



representation fj3.3[) the outer expansion. 

Let us substitute the outer expansion and (|3.2p into the eigenvalue equation 

H £ {k)4 j) =4 j) 4 j \ (3-5) 

and then group the coefficients at the like powers of - — . Considering the first two terms of the 

me 

(i) 
series obtained gives the following problems for the coefficients 4>\ . 



(( 






^+^-^T)K * ^ (-) 



^1 n tAO Aj) n -ph 

-^— = on T N ' , 4>\ J > = on T D , 

and the periodic boundary conditions (|1.4cp for the functions cpf are imposed at the side 
boundary of the domain H h . The conditions obtained split into two independent problems on 
Yr; and Il_. Their solvability in W^lL):) is equivalent to the orthogonality of the right-hand 
sides to the the functions (fr in ^(IIj.). It is easy to see that the solvability conditions will 

(i) (i) 

be satisfied only if a+ = or a_ = 0. On the other hand, in that case one cannot match the 

asymptotic expansions near the window ojq. Such difficulties are quite standard for problems 

involving a change of boundary conditions on a small part of the boundary, see e.g. [2Q|, Chs. 

2 and 6] or |251 Ch. M, §1], and they are usually overcome by considering the problem (|3.6[) in 

(i) 
a larger sense. Assume that <j)^ has a logarithmic singularity at the contraction point of the 

window ujq-. 

(j)[ j) (x) = b ( i ) In \x\ +Z < £ > + 0{\x\), x^O, ±x 2 >0, (3.7) 

(i) ~(i) (i) 

where b± , b± are some constants, and we assume that that the function cfr^ is sufficiently 

regular outside the point x = and satisfies all the boundary conditions mentioned above. 

Let us multiply the equations (|3.6p by <p and then integrate by parts in 11^. taking into 

account the asymptotics (|3.7p . This gives necessary solvability conditions for the problems 

(1Mb: 



nb if) = a U) f hfJt (j) d _ ^(27rm + ko)) 



(3.8) 



Using the standard arguments, see e.g. [25J Ch. M, §2, Th. 2.2] it is possible to show that these 

~(j) 
conditions are also sufficient. Note that there are no restrictions on the coefficients b± , and one 

can choose them freely by adding to cpf the functions of the form c±<p . 

The first terms of the expansions (|3.3p constructed above may serve as an approximation of 
the true eigenfunctions outside a certain neighborhood of the window Wq only, while the boundary 
conditions near the window are not satisfied. To construct the asymptotics near the window let 
us introduce the stretched coordinates £ = (Ci,^) = xe~ x and look for the asymptotics of the 
eigenfunctions in the form of an inner expansion: 

V, CO (x) = e 1 ^* 1 («, $» (0 + JUP* (0 + . . .) (3.9) 
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Let us substitute this ansatz and (|3.2p into the equation (|3.5[) and group the coefficients in front 

of e~ 2 and s~ 2 In - e in this formally obtained equality. This results in the following conditions 

(7) 
on the functions w p , p = 1,2: 

A,u#) = in IR 2 \r„ ^ = on F*, T* := O& \ [-1, 1]. (3-10) 

(i) 
To determine uniquely the functions w p we need to prescribe their behavior at infinity. This 

can be found from the matching of the inner expansion (|3.9|) with the outer expansion (|3.3|) . 

(i) (i) 

Let us represent the functions 0q and (^ in the variables £ in a neighborhood of 2 = 0: 

4> l i\x) + ^-^\x) = a < i ) +b { l ) 



lne 

2 } ln|£|+62 j ) + 0(e(|£| + l)), x^O, ±x 2 > 0. 



_^W)lv, ICI _L fcO')> 



lne 
By matching these relations with (|3.9p . we obtain 

w { j\0 = a { i ) + bf + o{l), |£|-+oo, ±6>0, (3.11) 

4 i) (0 = &± ) ln|£|+62' ) , lel^oo, ±&>0. (3.12) 

Therefore, one needs to find the solutions of the problems (|3.10p behaving at infinity in the 
prescribed way. 

The solvability of the problems (|3.10p is rather easy to study with the help of the methods 
of complex analysis. Namely, the domain K 2 \ T* can be conformally mapped onto the disk with 
the help of an explicitly given conformal map. The problem in the disk can be then studied in 
a standard way, one just needs to take into account the behavior at the image of infinity. Such 
an approach gives the following results on the solvability of the problems (|3.10p , (|3.1ip , (|3.12p . 

The problem (|3.10p with the infinity behavior (|3.1ip has the unique solution which is the 
constant one: Wq (£) = Co. This is only possible if the solvability conditions 

a (f) + b (J) = a if) +b (J) = CQ (3.13) 

hold. The problem (|3.10p with the condition (|3.12p can also be solved uniquely (up to a constant) 
under the compatibility conditions 

bf = -b { l\ (3.14) 

Its solution can be constructed explicitly: 

w[ j) (0 = b { l ) X(0 + c l , Cl = const, 

(3.15) 

X(£):=Reln(z + Vz 2 -l), z = £i + i&. 

Here the branches of the square root and the logarithm are chosen by the conditions \/T = 1> 
In 1 = 0. 

Let us emphasize some properties of the function X, which can be easily checked using 
its explicit definition. The function X is infinitely smooth in R 2 \ T* and is contiuous up to 
the boundary if this domain. It also has one-side derivatives on the set r*, which should be 
understood as a two-side cut. The function X is odd with respect to £2 and has the following 
asymptotic behavior at infinity: X(£) = ±(ln |£| +ln2) + 0(|£|" 2 ), |£| -)■ 00, ±£ 2 > 0. This last 
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equality leads to the condition f|3. 14|) . Due to the asymptotics (|3.12p . this equality also allows 
one to determine the constants b± : b± = b± In 2. 

By elementary computations one can easily check that the system (|3.13|) , (|3. 14[) is equivalent 
to the following one: 

_26«) - a { i ] + o^ = 0, 



n ,&') _ lr„U) . 07,0') 

d J+ ~d 



la^-^(a^+2b^) = 0. 



Let us express the numbers b± in terms of //0) using the formulas (|3.8p and substitute the 
expressions obtained into the last system of equations: 

(M - 2hn {j) )a {j) = 0, 

where 

2vrn + A; 27rm + & vr 

ft:=—^— , fc.= , C = 2 - (3-17) 

As we assume that the vectors a^> are non-zero, the quantities 2hp> 3 ' must be the eigenvalues 
of the matrix M, therefore, 

M CJ) (t) = £M f (3.18) 



/i(t, C) := m - /5 2 ) - (C + 1) + V(*(A + /3 2 ) + C - I)' + 4C, 



/ 2 (t, C ) := t^ - /3 2 ) - (C + 1) - ^{t(f3 1 + (3 2 ) + C-l) 2 +4C- 
The respective eigenvectors can be chosed by: 

^=(2^+1-^)' j = 1 > 2 ' (3 - 19) 

(i) 
The expressions obtained for a± completely determine the inner and the outer expansions for 

the eigenfunctions tpe • The constant c\ in (|3.15j) will be set equal to zero. Therefore, the formal 
construction of the asymptotics is completed. 

3.2 Justifying the formal asymptotics 

The first step of the justfication consists in a proof of the fact that the asymptotics (|3.2|) . 
(|3.3p deliver formal asymptotic solutions of the equation (|3.5|) . i.e. satisfy it with some small 

error. In our case one needs to obtain an error of order 0\ -, ttt ) . For that, it is necessary 

\\me\ z J 

to construct some additional terms in the outer and the inner expansions of the eigenfunctions 
ips and some additional summands in the asymptotics (|3.2j) . This is a standard situation in 
problems with singular perturbations. In our case, such a construction does not involve any 
conceptual difficulties but is rather technical. Furthermore, it consists essentially in a repeating 
of the similar constructions of [20], |25t Ch. H, §1]. Hence, we will not enter into details of the 
construction of additional terms, but prefer to formulate directly the final result. Let us denote 



WM =P i3 TT*i 



wK*) = e ' 



$\z) + ^/f\*))x(f! 
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(j) r \ /. f\x\ 

' ' IneJ \~ ^\e< 



E®(t):=Eo + 



+ K ; (0 + — ) (i-x 



where x = x( s ) is a smooth function which a equal to zero for s < 1 and equal to one for s > 2, 
and a is some constant. 

Lemma 3.1. The function ip £ ' belongs to C°°(H h ) n W 2 l (n h )- One /ias i/ie convergence 

$G) _). e^* 1 ^ in L 2 {U h ), E® -»■ £ - (3.20) 

There exists a number a € (0,1), functions ip£ (x,e) £ C°° (U h ) (IW2 (H h ) , and numbers E? (e) 
such that the functions 

V$ (x) := ®) (x) + # (x, e), £ £ ,* := ^ j) + M i} (e) 

satisfy the equation 

(H%k)-E^)$}=g¥\ (3.21) 

one /ias, uniformly in e and t, the estimates 

IWPMWwtmk) = 0(ln~ 2 e), \E®(e)\ = 0(ln" 2 e), (3.22) 



1 2 



ll#llL 2 (n^)=0(ln" 2 e). (3.23) 

Now let us consider the remainders of the asymptotics. By (|3.21|) and [261 Ch. V, §3.5, Eq. 
(3.21)] there holds 

2 / (3) ,i,(p)\ 
p =i E e — E £ 

where, E £ ' are the eigenvalues of the operator H e (k) converging to Eq as e — > +0, ?/4 are 
the respective eigenfunctions which are supposed to be orthonormal in ^(n ), and R £ is the 
reduced resolvent considered as an operator in L2(H h ) and satisfying the estimate 

lim = d^(E^\a(Hs { k e ))\{EP,EP}) * C ' 

(i) 
where the constant C does not depend on e and t. Furthermore, the operator R £ acts in the 

orthogonal complement of {ip £ ,ip £ } in £2(11 ). 
Let us denote 

tJj) j,(p)\ 

. , , \9e ,Ve L 2 (Il h ) . , ■. , . , ■. . , ... 

E £ — E e 

It follows from the properties of the operators R £ described above and from (|3.24p . (|3.23jl . 
d32QD, (I3TT9J1 . lETT7|) . (IXTHjl . (ETHl that 

II^SllLcn*) = E l^l a + ll4 j) #llL ( n^ i = I, 2, 
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(^ 1) ,^ 2) )L 2 (n^)=E^P^ + (4 1) 5«,4 2) #k 2 (n^) J 
P =i 

and 

IIA/lk 3 = ll^lkcnfc) + o(l), e -> +0, 

(Ai,A 2 ) R2 = (^ol^L.m +o(l) = o(l), e^ +0, 

uniformly in t £ [—to, to]. At cost of redenoting ?/4 by ?/4 and vice versa, due to the above 
relations, we conclude that one has the estimates 

|Aii|>C>0, |,4 22 |^C>0, 

which are uniform in e and t G [—to, to]- Therefore, by the definition of the functions Aj p and 
by the estimate (|3.23|) , we obtain the estimates of the remainder 

E®(t)-E®{t) = OQn- 2 e) 

uniformly in t E [—to, to]- Now, taking into account (|3.22|) . we obtain the final form of the 
asymptotics for the eigenvalues Eg 1 uniformly in t € [—to, to]: 

E^(t) = E + ^^ + O(ln- 2 e), j = 1,2. (3.25) 



At the end of the section let us remark that, using the methods of the works cited above, 
one can construct an asymptotics of the eigenvalues Eg up to a certain power of e by summing 
up the logarithmic part. On the other hand, such a precision is redundant for our objectives. 



3.3 Studying the asymptotic formulas 

Let us proceed with the qualitative analysis of the band functions near k = k^ using the asymp- 
totics obtained. To start, let us determine the extrema of the functions fj in (|3. 18[) . 

It follows from the inequality (|1.7a|) that the numbers f3\ and /3 2 in (|3.17|) have the same 
sign. Hence 

|0i - H = ||£i| - N, |/3i +/3 2 | = |/3i| + |&|, 

which implies |/3i — /3 2 | < \fi\ + /3 2 |- Therefore, 

lim /i(t,C) = +oo, lim f 2 (t, Q = -oo. 

t— >±oo t— >±oo 

Using the above and the continuity of the functions f\ and / 2 with respect to t one sees that 
they attain their respective global minimum and maximum at some finite points. Calculating 
these points in the standard way we obtain 

m m/l(i,C) = /l(imin,C), max/ 2 (t,C) = /2(*max,C)> 



2/3 / C 1-C 



Px + fa V 1 - /3 2 Pi + P. 



2 



2/3 / C + 1-C 



/3i + /3 2 V 1 " /3 2 /3i + /3 2 ' 
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where 

a .__ 01 ~^ __ 7r(m + n) + k 
Pi + @2 n(n - m) 

The values of the global minimum of the function f± and of the global maximum of the function 
J2 are given by the formula 



/i(W, C) = 2VCU-/3 2 ) - /3(C - 1) - C - 1 = r+, 
/ 2 (W, = -2 X /C(1-/? 2 ) - /3(C - 1) - C - 1 = r_. 

The inequality lne < together with the preceding relations and the equalities <|3.18|> . (J3.25D . 
(1331) give 

max £?W =E - —rr. , min /i(t, C) + of 



t6[-to,*o] 4/i| lne| te[-t ,t ] V|lne| 2 



4/i|lne| \|lne| 2 

min Ej® =E - — - f max / 2 (t, Q + o( 



te[-t ,t ] 4h\ In e\ t&[-t , t ] ' \|lne| 

1 



~-Eo ~ -TTT^ + O 



4/i|lne| \|lne| 

where the parameter to can be chosen sufficiently large to have t m \ n £ [— to, to], £ max £ [— to, to]. 
It follows from the last equality and from ()3.26|) that, for e sufficiently small, 

min Ef\t)- max #(t)= r )"^ +o(— — 
te[-to,to] te[-to,*o] |lne| V|lne 



/i|lne| V|lne| 2 7 

Taking into account the inequality (|3.ip and choosing the parameter A sufficiently large we 
conclude that there is a gap in the spectrum of the operator H £ , 

Eo- TT ^ + o( TT ^),Eo-^^ + o( l 



4ft,|lne| \|lne| 2 /' 4/i|lne| \|lne| 2 

which proves (|1.8|) and (|1.9p . It remains to show the validity of the asymptotics (jl.lOp . 

As the function f\{t, Q) attains its global maximum at the point t m j n only, for all t at a finite 
distance from t m \ n there holds f\(t, Q) > /i(£ m in, C)- Therefore, the point of the maximum of the 
function E[ £ should tend to t m ; n as e — > +0. As 

fl(t, C) = /l(imin, C) +0{(t- t min ) 2 ), t -> t min , 

the equality 

^ ll(4) = E »-w^ + °(n^) 

can only be satisfied for t situated at the distance of at most O ( | In s\ ' 2 ) from i m i n . This gives 
the asymptotics ([l.iup for kt(e). The asymptotics for k^(e) is proved in the same way. 
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